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1.0 Configuration Determination 


To determine the best configuration from all candidate configurations, it was necessary 
first to calculate minimum system weights and performance To optimize the design, 
necesary to vaiy configuration-specific variables such as total system weight, thrust t 
SSS5 durations, total thrust available, and mass fraction for the system. 
Optimizing each of these variables at the same time is technically unfeasible and not 
necessarily mathematically possible. However, discrete sets 

which will eliminate many candidate configurations. From the most promising 
remaining designs, a final configuration can be selected. 

To ease this optimization process an interactive spreadsheet 

thrust to weight ratios, number of main engines, vehicle diy weight , andjnitial (p e 
separation) burn time. From this information, it calculated maximum weights 
propellant mass fractions available for each vehicle and the final post-separation bum 

time for the orbiter. 

After selecting the number of engines and thrust to weight ratios for both the booster and the 
orbiter the length of the initial bum was varied to maximize fuel available to the orbiter for 
the second phasebum. This maximum weight was then used to determine tank sizings 
which dictated total volume necessary within the booster and orbiter. 

Included are the three meet important designs considered, oh . which eta* '■W™^*** 
the design criteria set forth in a Marshall Space Flight Center study of the Shuttle II, the 
configuration used in the initial proposal; and the final configuration presented an 
analyzed in this report Note that these three systems represent important iterative steps in 

our final design. 

A listing by cell of the formulas used to generate the aforementioned data is included for 
reference. 
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ORBITER DAT* 


.Thrust /Weight 

1.2000 

Initial Boost Engines 

o 

o 

Total Thrust at Takeoff 

1, 454 , 100 . 00 lbf 

Available Takeoff Weight 

1, 211, 750 . 00 lbm 

Dry Weight 

180,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

20,000.00 lbm 

Total Orbiter Weight 

260, 000 . 00 lbm 

Maximum Fuel Weight 

951,750.00 lbm 

Fuel Available Initially 

951, 750 . 00 lbm 

Initial Boost Engines 

4 . 00 

Initial Fuel Consumed 

505, 080 . 80 lbm 

Primary Fuel Available 

446, 669.20 lbm 

Auxilliary Fuel 

496, 952 . 57 lbm 

Total Fuel Available 

94 3, 621.. 7 7 lbm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1,740,000.00 lbf 

Approximate Final T/W 

1.4456 

BOOSTER DATA 


Thrust /Weight 

1.4000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1, 817, 625.00 lbf 

Available Takeoff Weight 

1,298,303.57 lbm 

Dry Weight 

160,000.00 lbm 

RCS Fuel (MMH/N04) 

10,000.00 lbm 

, Total Booster Weight 

170, 000 . 00 lbm 

Fuel Weight 

1, 128, 303 . 57 lbm 

Fuel Available Initially 

1,128,303.57 lbm 

Initial Boost Engines 

5.00 

Initial Boost Fuel 

631,351.00 lbm 

Reservoir Fuel Pumped 

496, 952 . 57 lbm 

OVERALL DATA ------- 


Thrust/Weight 

1.3034 

Initial Boost Engines 

9 . 00 

Total Thrust at Takeoff 

3, 271,725 . 00 lbf 

Available Takeoff Weight 

2, 510, 053 . 57 lbm 

Dry Weight 

340,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

30,000.00 lbm 

Total Dry Weight 

430, 000 . 00 lbm 

Fuel Weight 

2, 080, 053 . 57 lbm 

Fuel Available Initially 

2,080, 053 . 57 lbm 

Initial Boost Engines 

v> 

o 

o 

Initial Boost Fuel 

1, 136, 431 . 80 lbm 

Fuel Available Finally 

943, 621.77 lbm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


133.00 sec 

248.48 sec 

381.48 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 Ibm/sec 
0.0700 ft~3 
1.1490 ft *3 
62.43 lbm/ f t A 3 

135, 964.29 lbm 
815, 785.71 lbm 
31, 112 . 40 ft *3 
11, 372 . 68 ft ^3 
42 , 485.07 ft ~3 
88.0301 ft 
21.4203 ft 
109 .4504 ft 


161,186.22 lbm 
967,117.35 lbm 
36,883.88 ft "3 
13,482.36 ft ^3 
50,366.23 ft * 3 
69.4704 ft 
25.3939 ft 
94.8643 ft 
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ORBITER DATA 


Thrust /Weight 

1.2000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1, 817 , 625 . 00 lbf 

Available Takeoff Weight 

1, 514 , 687 . 50 lbm 

Dry Weight 

180 , 000.00 lbm 

Payload 

60, 000.00 lbm 

RCS Fuel (MMH/N04) 

20,000.00 lbm 

Total Orbiter Weight 

260 , 000.00 lbm 

Maximum Fuel Weight 

1, 254 , 687 . 50 lbm 

Fuel Available Initially 

1, 254, 687 .50 lbm 

Initial Boost Engines 

5 . 00 

Initial Fuel Consumed 

564,893.00 lbm 

Primary Fuel Available 

68 9, 794 . 50 lbm 

Auxilliary Fuel 

563, 410 . 57 lbm 

Total Fuel Available 

1, 253 ,205.07 lbm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1,740,000.00 lbf 

Approximate Final T/W 

1.14 99 

BOOSTER DATA 


Thrust/Weight 

1.4000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1, 817, 625 . 00 lbf 

Available Takeoff Weight 

1, 298, 303 .57 lbm 

Dry Weight 

160,000.00 lbm 

RCS Fuel (MMH/N04 ) 

10, 000.00 lbm 

Total Booster Weight 

170, 000 . 00 lbm 

Fuel Weight 

1,128,303.57 lbm 

Fuel Available Initially 

1, 128, 303 . 57 lbm 

Initial Boost Engines 

5.00 

Initial Boost Fuel 

564,893.00 lbm 

Reservoir Fuel Pumped 

563, 410 . 57 lbm 

OVERALL DATA 


Thrust/Weight 

1.2923 

Initial Boost Engines 

o 

o 

o 

H 

Total Thrust at Takeoff 

3, 635, 250 . 00 lbf 

Available Takeoff Weight 

2 , 812 ,991.07 lbm 

Dry Weight 

340,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04 ) 

30 , 000.00 lbm 

Total Dry Weight 

430 , 000.00 lbm 

Fuel Weight 

2 , 382 ,991.07 lbm 

Fuel Available Initially 

2 , 382 ,991.07 lbm 

Initial Boost Engines 

10 . 00 

Initial Boost Fuel 

1,129,786.00 lbm 

Fuel Available Finally 

1, 253 ,205.07 lbm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq* Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


119.00 sec 

330.00 sec 

449.00 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 ibm/sec 
0.0700 ft A 3 
1.1490 ft A 3 
62.43 lbm/ft~3 

179, 241.07 lbm 
1, 075 , 446.43 lbm 
41,015.32 ft * 3 
14,992.55 ft * 3 
56,007.87 ft ^3 
116.0496 ft 
28.2383 ft 
144.2879 ft 


161, 186.22 lbm 
967,117.35 lbm 
36,883.88 ft *3 
13,482.36 ft *3 
50,366.23 ft *3 
69.4704 ft 
25.3939 ft 
94.8643 ft 
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ORBITER DATA 


Thrust/Weight 

1.2000 

Initial Boost Engines 

4.00 

Total Thrust at Takeoff 

1, 454,100.00 lbf 

Available Takeoff Weight 

1, 211,750 . 00 lbm 

Dry Weight 

180,000.00 lbm 

Payload 

60, 000 . 00 lbm 

RCS Fuel (MMH/N04) 

20,000.00 lbm 

Total Orbiter Weight 

260,000.00 lbm 

Maximum Fuel Weight 

951, 750.00 lbm 

Fuel Available Initially 

951,750.00 lbm 

Initial Boost Engines 

4 . 00 

Initial Fuel Consumed 

562. 044.80 lbm 

Primary Fuel Available 

389, 705 .20 lbm 

Auxilliary Fuel 

544, 897 . 09 lbm 

Total Fuel Available 

934, 602.29 lbm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1,740,000 . 00 lbf 

Approximate Final T/W 

1.4566 

BOOSTER DATA 


Thrust/Weight 

1.4000 

Initial Boost Engines 

6.00 

Total Thrust at Takeoff 

2, 181,150.00 lbf 

Available Takeoff Weight 

1, 557, 964.29 lbm 

Dry Weight 

160,000.00 lbm 

RCS Fuel (MMH/N04) 

10,000.00 lbm 

Total Booster Weight 

170, 000 . 00 lbm 

Fuel Weight 

1,387,964.29 lbm 

Fuel Available Initially 

1,387, 964.29 lbm 

Initial Boost Engines 

6.00 

Initial Boost Fuel 

843, 067 . 20 lbm 

Reservoir Fuel Pumped 

544,897.09 lbm 

OVERALL DATA 


Thrust/Weight 

1.3125 

Initial Boost Engines 

10.00 

Total Thrust at Takeoff 

3, 635,250 . 00 lbf 

Available Takeoff Weight 

2, 769,714 .29 lbm 

Dry Weight 

340, 000. 00 lbm 

Payload 

60, 000 . 00 lbm 

RCS Fuel (MMH/N04) 

30,000.00 lbm 

Total Dry Weight 

430, 000 . 00 lbm 

Fuel Weight 

2, 339, 714 .29 lbm 

Fuel Available Initially 

2,339,714 .29 lbm 

Initial Boost Engines 

10 . 00 

Initial Boost Fuel 

1,405,112.00 lbm 

Fuel Available Finally 

93 4,602.29 lbm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


148.00 sec 

246.10 sec 

394.10 sec 

363 , 525.00 Ibf 
435 , 000.00 lbf 

949.40 lbm/sec 
0.0700 ft *3 
1.1490 ft ^3 
62.43 lbm/ft~3 

13 5, 964 . 29 lbm 
815,785.71 lbm 
31, 112 .40 ft ^3 
11 , 372 . 68 ft ^3 
42,485.07 ft*3 
88.0301 ft 
21 . 4203 ft 
109.4504 ft 


198,280.61 lbm 
1, 189, 683 . 67 lbm 

45.372.10 ft A 3 

16.585.10 ft~3 
61, 957 .20 ft ~ 3 

85.4579 ft 
31 . 2379 ft 
116.6958 ft 
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ORBITER DATA 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
Payload 

RCS Fuel (MMH/N04) 

Total Orbiter Weight 
Maximum Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Fuel Consumed 
Primary Fuel Available 
Auxilliary Fuel 
Total Fuel Available 
Final Boost Engines 
Approximate Final Thrust 
Approximate Final T/W 

BOOSTER DATA 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
RCS Fuel (MMH/N04) 

Total Booster Weight 
Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Boost Fuel 
Reservoir Fuel Pumped 

OVERALL DATA 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
Payload 

RCS Fuel (MMH/N04 ) 

Total Dry Weight 
Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Boost Fuel 
Fuel Available Finally 


1.2 

5 


-B3*E6 

lbf 

-B4/B2 

lbm 

180000 

lbm 

60000 

lbm 

20000 

lbm 

-B6+B7+B8 

lbm 

-B5-B9 

lbm 

-B10 

lbm 

-B3 


— B13*18*X2 

lbm 

-B12-B14 

lbm 

— B3 4 

lbm 

-B15+B16 

lbm 

4 


-B18*E7 

lbf 

— B1 9/ (B9+B17) 


1.4 


5 


— B23 *Z6 

lbf 

-B24/B22 

lbm 

160000 

lbm 

10000 

lbm 

-B2 6+B27 

lbm 

— B25-B28 

lbm 

— B2 9 

lbm 

-B23 


— B3 2 *18*12 

lbm 

-B31-B33 

lbm 

-(B4+B24) / (B5+B25 ) 


-B3+B23 


-B4+B24 

lbf 

-B38/B36 

lbm 

-B6+B26 

lbm 

1 

0> 

lbm 

-B8+B27 

lbm 

•B40+B41+B42 

lbm 

•B10+B29 

lbm 

-B12+B31 

lbm 

-B37 


■ B4 7 *E2 *E8 

lbm 

-B46-B48 

lbm 



Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 



Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


.82 9/7 Ibm 

— B2 9*6/7 lb® 

-127 / (19*111) ft'3 

-*28/(110*111) ft'3 

—12 9+130 ft~3 

—129/(910/4*28*2) ft 

-130/(910/4*28*2) ft 

— *32 +133 ft 
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2.0 STME Performance Calculation Summary 

A model of the total thrust available from the Space Transportation Main Etipnes used on 
rtw. Shuttle II was necessary to evaluate its launch profile. Thrust is a function of many 
^riabte^ hiclutoig^^^e^ ixpansion ratio, exit velocity, ambient pressure and exit 
pressure Since the STME uses a translatable skirt, the expansion ratio of the engine 
be set at 'discrete values of 35:1, 70:1, and 150:1. To insure that the engine * ope 
efficiently during all points in the launch, it was necessaiy to calculate the altitudes at 
which the skirt should be translated. This analysis was performed by calculating th 
coefficients for all three nozzles and finding the points of intersection. 

These calculations were based on the data presented in the Aeroject contractor's report for 
^ STstandard atmosphere was used as a model for pressure and dens.ty as 
functions of altitude. Most of the engine data was solved using simple, iterative 
techniques, isentropic relations, and graphical methods. 


2.1 Analysis 

Given an expansion ratio for a particular nozzle and an appropriate value for the 
coefficient of specific heats, the exit Mach number may be determined. 


Assuming y = 1.24 for a mixture of liquid oxygen and liquid hydrogen in a 6:1 mass ratio 
and the first (fixed) nozzle expansion ratio of 35:1, this isentropic relation holds: 

At _ A— j I +■ 1*1 e- 1- "1 j z I t'H 

/X* **«-L J 

This function may be solved iteratively to obtain exit Mach numbers for the nozzles. This 
was solved using the included iterative, and the following were obtained as the exit Mach 
numbers: 


Nnasrie Expansion Ratio 
35:1 
70:1 
150:1 


Frit Mach Number 
4.369 
4.931 
5.581 


Next, the flow in a nozzle may be assumed to be very nearly isentropic, such that the 
chamber pressure may be related to the nozzle exit pressure through the isentropic relation 

Given a chamber pressure of 2385 psia under normal power levels and 3125 psia under 
emergency conditions, the following exit pressures were obtained: 


Pe (NFL) PeflSPLl 
Fore = 35:1 5.069 6.642 

For e = 70:1 2.058 2.696 

Fore =150:1 0.771 1.010 


For a given expansion ratio, a nozzle will have a maximum efficiency at a certain design 
altitude. At this altitude, the exit pressure is equal to the ambient pressure of the 
atmosphere. From the 1962 standard atmosphere, the design altitudes were determined to be 

hd (Normal P ower Level) (EmerggnCY POWCT Lgrel al 

For e = 35:1 26,700 ft 20,400 ft 

For e = 70:1 45,800 ft 40,200 ft 

For e = 150:1 66,300 ft 60,600 ft 


In order to obtain the proper altitudes at which the sections of the nozzle should be extended, 
the following procedure is followed. 


The thrust coefficient is defined as: 

Or = 


Pt a* 


Since thrust varies with altutde, this expression may be rewritten as. 
-- 



w 


Pr 


2 


Ac 
A * 
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For example, consider the nozzle with an expansion ratio of 70:1 at an 


altitude of 50,000 feet. 


Cj ■■ * (j^f) 


11 355 


2-050 * 


- ( 10 ) 


CIS 


Similarly, thrust coefficients for other nozzle and atmospheric conditions may be 
calculated. As shown in Figure 2.1.1, the points where these graphs cross indicate the 
altitudes where the nozzles should be translated for maximum efficiency. 

For Normal Power Levels, the translation altitudes were 56 ’ 5 °° 

feet For Emergency Power Levels, the altitudes were 31,100 feet and 50,800 feet. 
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21000 


22000 


23000 


24000 


25000 


26000 


27000 


28000 


29000 


30000 


31000 


32000 


33000 


34000 


35000 


36000 


37000 


38000 


39000 


40000 


41000 


42000 


43000 


44000 


45000 


46000 


47000 


48000 


49000 


50000 


51000 


52000 


53000 


54000 


55000 


B 


P 


14.7000 


14 . 1700 


13 . 6600 


13.1700 


12.6900 


12.2300 


11.7800 


11.3400 


10.9200 


10.5000 


10.1100 


9.7200 


9.3460 


8.9840 


8.6330 


8.2940 


7.9650 


7.6470 


7.3390 


7.0410 


6.7540 



UP h , 1*3 5 


1.6185 


1.62631 


1.6338 


1.64101 


1.6480 


1.6548 


1.6614 


1.6679 


1.67401 


1.6802 


1.685 91 


6.20701 


5.94701 


5.6960 


5.45401 


5.2200 


4 . 994 ol 


4.7770 


4.5670 


4 .36401 


4.1690 


3,9810 


3,80001 


3.6260 


3.4580 


3.29701 


3.1420 


2 .99401 


2.8540 


2.7200 


2.5920 


2.4710 


2 .35501 


2.1390 


2 . 0 390 l 


1.9430 


1.8520 


1.7650 


1.6820 


1.60301 


1.52801 


1.4560 


1.3880 


1.3230 


1.6971 


1.7024 


1.70761 


1.71261 


1.7174 


1.7220 


1.72661 


1.7309 


1.7352 


1.7392 


1.7432 


1.74701 


1.7507 


1.75421 


1.7577 


1.7610 


1 .7642 


1.7672 


1.7702 


1.7731 


1.7758 


1 .77851 


1.7811 


1.7835 


1.78591 


1.7882 


1 .79031 


1.7924 


1.79431 


1.7962 


1 .79801 


1.7997 


1.8013 


1.8029 


1.80431 


1.8058 


1.8071 


1.8084 


1.80961 


1.8107 


1.8118 


1.8129 


1.8139 


1.81491 


1 . 6696 


1.6756 


1.6813 


1.6868 


1.6921 


1.6973 


1.7023 


1.7073 


1.7120 


1.7167 


1.7210 


1.7254 


1.7296 


1.7336 


1.7376 


1.7414 


1.7451 


1.7486 


1.7521 


1.7554 


1.7586 


1.7617 


1.7647 


1.7677 


1.7705 


1.7732 


1.7758 


1.7783 


1.7808 


1.7831 


1.7854 


1.7876 


1.7897 


1.7917 


1.7937 


1.7955 


1.7973 


1.7991 


1.8007 


1.8023 


1.8038 


1.8052 


1.8066 


1.8079 


1.8091 


1.8103 


1.8114 


1.8125 


1.8135 


1.8145 


1.8154 


1.8163 


1 . 8172 


1.8180 


1.8187 


1.8194 


1.4649 


1.4798 


1.4942 


1.5083 


1.5218 


1.5350 



1.5840 


1.59551 


1.6064 


1.6171 


1.6274 


1.63731 


1.6470 


1.6563 


1.66531 


1.68251 


1.69071 


1.6986 


1.7062 


1.7136 


1 .72071 


1 .73421 


1.7405 


1.7467 


1.75271 


1.7584 


1.7639 


1 .76921 


1.7793 


1.7840 


1.78851 


1.7929 


1.79701 


1 .80091 


1.8047 


1.8082 


1.81161 


1 .81491 


1.8180 


1 .82091 


1.8237 


1.8264 


1.8289 


1.8314 


1.8337 


1.83591 


1.8380 


1.8400 


1 .84191 


1.5515 


1.5633 


1.5748 


1.5857 


1 . 5965 L 


1 .60681 


H 


£■150 


1.2164 


1.6267 


1.6361 


1.6455 


1.6543 


1.66301 


1.6795 


1.6874 


1.6950 


1.7023 


1.7094 


1.7163 


1 . 7295 ) 


1.73571 


1.7532 


1.7586 


1.7638 


1.76891 


1.7784 


1.7830 


1.7874 


1.7916 


1.7956 


1.7995 


1.8033 


1.8069 


1.8104 


1.8137 


1.8168 


1.8198 


1.8227 


1.8254 


1.8280 


1.8305 


1.8328 


1.8351 


1.8372 


1.8393 


1.8412 


1.8431 


1.8448 


1.8465 


1.8481 


1.8497 


1.8511 


1.0308 


1.0629 


1 .0937 


1.1239 


1.1529 


1.1812 


1.2088 


1.2352 


1.2617 


1.2862 


1.3107 


1.3342 


1.3570 


1.3791 


1.4004 


1 .4605 


1.4792 


1 .4973 


1 .5148 


1.5317 


1.54 80 


1.5638 


1.5790 


1.5937 


1.6079 


1.6216 


1. 6348 


1.6476 


1 . 6598 


1.6717 


1. 6830 


1.6940 


1.7046 


1.75101 


1.7590 


1.7666 


1 . 7809 ) 


1 .7875 


1 . 7938 


1 .7998 


1.8056 


1.8110 


1.8163 


1.8212 


1 . 8259 


1 . 8305 ! 


1.8347 


1 . 8388 


1.2899 


1.3129 


1.3350 


1.3566 


1 .397 9 


1.4 180 


1 .4 368 


1.4555 


1.4 908 


1.5077 



1 .5550 


1.5698 


1.584 1 


1.5978 


1.6241 


1 . 6366 


1.6486 


1 . 6602 


1.6823 


1 .6927 


1 .702 8 


1.7126 


1 .7 309 


1.7396 


1 .7 480 



1 .7783 


1 .7850 


1 .7915 


1.7976 


1 . 8034 


1 .8090 


1.8143 


1 .8194 


1.8242 


1 . 8288 


1.8331 


1.8373 


1.8413 
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Figure 2.1 .1 : Plot of Thust Coefficients vs. Altitude (NPL) 
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Figure 2.1.2: Thust Coefficients vs. Altitude (EPL) 
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3.0 Mission Analysis 

To prove the flight worthiness and technical feasibility of the ASTS as proposed, a launch 
profile was constructed and simulated. This simulation consisited of a dynamic model of 
the vehicle’s instantaneous position, velocity, acceleration, thrust, weight, orientation, 
and flight path with respect to the horizon. In addition to this vehicle-specific data, 
atmospheric conditions such as temperature, pressure, and density were modeled. 

The basic format of this simulation was a fourth-order Runge-Kutta integration scheme. A 
number of asumptions were necessary to develop this model of the launch trajectory. First, 
the model is constrained by a run time of approximate seven minutes and a total distance 
downrange not exceeding 800 nautical miles, so the surface of the earth was modeled as 
being flat. In addition, rotational effects of the earth were neglected. This assumption 
cannot degrade the validity of results, since it could only create conditions more favorable 
to a due east launch from Kennedy Space Center. Drag forces were neglected for two 
principal reasons: they can be considered small in comparison to thrust forces necessary 
for a ballistic trajectory, and they would not exceed the advantages gained by 
incorporating the earth's rotation. 

All initial conditions were set based on values calculated from the design process outlined 
in Section 1.0. Bum times, initial weights, total thrust available, and fuel available are the 
the most noteworthy examples. Thrust was computed as a funciton of altitude, based on 
translation altitudes for the extendable skirt on the STME nozzles, pressure forces, and 
known values of sea level and vacuum values. It was estimated that the maximum 
acceptable pitch rate was 0.5 degrees per second to retain stability during launch. This is 
the value used in the program for all pitch manuevers. 

Before finalizing the launch profile, it was necessary to establish the range of orbits into 
which the orbiter could be successfully launched. Based on current STS mission analyses 
and estimations of maximum altitude and velocity available at this altitude, our launch 
goal was projected to be an altitude of 100 nautical miles and a velocity equal to that of the 
local circular speed. Inclination was varied during flight both to achieve this orbit and to 
place the booster in an alignment practical for separation. 


program AE447Prog; 

($MC68020+} 

{SMC68881+} 

{SSETC Elems881:-true) 


Written by Del Johnson (in TML Pascal II vl.00 under MPW v2 .02) 

Written”for AE447-AE449 Shuttle II design group (Fall 19 88 - Spri ng 198in 


This program is designed to model the flight of design configurations for the 
Shuttle II Launch System. Design data was used in the model- *® S f. tune the 
execution of these programs was used in an iterative process to fine tune the 

model as presented in the report . 


Definition of variables u[i] : 


( Integral variables- } 


u [ 1] 
u( 2] 
u[ 3] 
u [ 4] 


u[ 
u ( 


5] 

6 ] 


- u[x] - Location along the x-axis 
» u (xd] — Time rate of change of location along the 
m u [y] — Location along the y-axis 

„ u [yd] ■ Time rate of change of location along the 
m u[th] — Angle with respect to the y-axis 
« u[thd] « Time rate of change of angle with respect 


x-axis 

y-axis 

to the y-axis 


( Other variables— } 
u ( ] - u [alpha]- Angle of attack 

u[ ] - u [gamma]- Flight path angle 

u [ ] - u[wt] - Total vehicle weight 

u [ ] - u[wtd] - Time rate of change of total vehicle weight 

u ( ] - u[dens] - Density of the atmosphere 

u[ ] - u [press] - Density of the atmosphere 
u [ ] - u[ thrust] - Thrust produced by all engines firing 


*) 

uses MemTypes , SANE ; 


const 


{ Number of equations- } 

numlnt Equations - 6; 

numTotEquations - 11; 


{ Time Considerations- ) 

tbegin - 0.00; [ seconds } 

tend - 420.00; { seconds } 

dt - 0.10; { seconds ] 

cint - 10? 

cols - 4; 

{ Subscripts of the variables of motion— } 


) 
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xd “ 1 * 

x ■ 2; 

yd “ 3 ; 

y “ 4; 

wt “ 5 ; 

th - 6; 

vtot ” 7 ; 

gamma “ 8 ; 

thrust “ 9 ; 

press ” 10; 

dens ” 11; 

Separation “ 1 
Shutdown ” 2 


StartPitchl - 3; 

StopPitchl ■ 4; 

StartPitch2 - 5; 

StopPitch2 • 6 ; 

{ Constants related to the problem.- ) 

{ Mathematical constants.- } 

conv - 57.29577951; { degrees /radians } 


( 


Physical 

Re 

gO 

rho 

rhoO 

Psl 


constants.- } 

- 20938912.00000; ( feet } 

- 32.17400; { feet/second~2 } 

- 0.0023768; { pound-mass /feet *3 ) 

- 0.00237 68; { pound-mass/feet^ } 

- 14.70000; { pound-force/ inches / '2 l 


{ Vehicle characteristics- ) 

thrO - 363535.00000; { lbf } 

wt d0 - -949.40000; { lbm/sec } 


pressAlt “ 25000 
dens Alt “ 25000 

Ael - 3247 
Ae2 - 6503 
Ae3 - 14526 


00000; { feet } 
00000; { feet } 

22085; { inchest ) 
88219; { inches A 2 } 
72443; { inchest } 


designAlt 1NPL 
designAlt2NPL 
designAlt 3NPL 
designAlt 1EPL 
designAlt 2EPL 
de s ignAl t 3 EP L 


26700.00000; { feet } 
45800.00000; { feet } 
66300.00000; { feet } 
20400.00000; { feet } 
40200.00000; { feet } 
60600.00000; { feet } 


transAltlNPL - 36900.0000; { feet } 
transAlt2NPL - 56500.0000; { feet } 
transAltlEPL - 31100.0000; { feet } 
transAlt2EPL - 50800.0000; ( feet } 


PelNPL - 5.069; { lbf/feet~2 } 
Pe2NPL - 2.058; { lbf/feet~2 } 
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Pe3NPL - 0.771; { lbf/feet A 2 } 
PelEPL - 6.642; { lbf/feet A 2 } 
Pe2EPL - 6.296; { lbf/feet A 2 } 
Pe3EPL - 1.010; { lbf/feet A 2 } 


thrMaxlNPL * 411259.15000; 
thrMax2NPL - 424577.04000; 
thrMax3NPL - 435000.00000; 

thrMaxlEPL * 545789.15000; 
thrMax2EPL - 556312.00000; 
thrMax3EPL - 580000.00000; 

{ Problem intial conditions... } 


xdO 

- 

0.00 

{ 

feet /second 

1 

xO 

• 

0.00 

{ 

feet } 


ydO 

- 

0.00 

{ 

feet/second 

} 

yO 

- 

0.00 

{ 

feet } 


thO 


90.00 

{ 

degrees } 


vtotO 

- 

0.00 

{ 

feet /second 

} 


type 

uArray 
name Array 
print range 


array 
array 
1. .32; 


[1. .numTotEquations] 
[1. .numTotEquations] 


of extended; 
of str255; 


var 


outfile : text; 
tab : char; 

3 p : array [1 . .numTotEquations] of integer; 

DryWtBooster, DryWtOrbiter, wtO : extended; 
t : extended; 

time : array [1..6] of extended; 

Jc ; array [ 1 . . 4 , 1 . .numlntEquations] of extended; 

u,uu : uArray; 

i, j,m,n : integer; 
printset : set of print range; 


{ Variables of motion^ } 

wtd, thr, thd : extended; 

thrMax : extended; 


function g(alt : extended): extended; 
begin 

g g0*sqr (Re) /sqr (Re+alt ) ; 

end; 


function density (alt : extended): extended; 
begin 

density :* rho0*exp {-alt /dens Alt) ; 

end; 

function pressure (alt : extended): extended; 
begin 

pressure :« Psl*exp (-alt/pressAlt ) ; 

end; 
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function pitchRate(t 
begin 


extended) : extended; 


(t < Time [StartPitchl] ) then pitchRate :« 
if (t < Time [StopPitchl] ) then pitchRate 

if (t < Time [StartPitch2] ) then pitchRate 

if (t < Time [StopPitch2] ) then pitchRate 

if (t < tend) then pitchRate 0.00; 


end; 

function thrustAvail (alt 


extended) : extended; 


0.00 else 

-0.50 else 
:« 0.00 else 

-0.50 else 


var 

Ae,thrl : extended; 

begin 

if (alt<transAltlNPL) then begin 
Ae Ael; 
thrl :» thrMaxlNPL; 


end 

else 

if (alt<transAlt2NPL) then begin 
Ae :■ Ae2; 
thrl :» thrMax2NPL; 

end 

else begin 

Ae :• Ae3; 

thrl :« thrMax3NPL; 

end; 

thrustAvail :* thrl - pressure (alt) *Ae; 

end; 

procedure initVariables ; 
begin 

t : - tbegin ; 
u[yd] :» ydO; 
u [y] yO; 

u[xd] xdO; 
u[x] xO; 

u[wt] :* wtO; 
u[th] :» thO; 
u [gamma] :» 90; 

u[vtot] 2 - sqrt (sqr (u [yd] ) +sqr (u[xd] ) ) ; 
u [press] pressure (u [y] ) ; 

u[dens] density (u [y] ) ; 

end; 

procedure print (style : integer); 
var 

i : integer; 
begin 

case style of 
1 : begin 

write (outfile, f t - v ,t:6:2); 

for i :■ 1 to numTotEquations do begin 

if (i in printset) then write (outfile, 1 u[ v ,i:2, f ] 

- 1 ,u [i] : 1 2 : 2 ) ; 

if (i in printset) then write ( 1 u[ f ,i:2, f ] 

- f # u[i] :12:2); 

end; 


} 
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') ; 
') ; 


writeln {outf ile) ; 
writeln; 
end; 

2 ; begin 

i 0; 

writeln (outf ile, 'Time - ' ,t:7:2); 
writeln ( 'Time m ’,t:7:2); 

for j :■ 1 to numTotEquations do 
if (j in print set) then begin 
i i+1; 

write (outf ile, 'u[\ j:2, '] ,u[ j] : 12 :2, ' 

write ( 'u[',j:2,'] “',u[j]:12:2,' 

if (i mod cols - 0) then begin 
i :« 0; 

writeln (outf ile) ; 
writeln; 
end; 
end; 

if (i<>0) then begin 
writeln (outf ile) ; 
writeln; 
end; 
end; 

3 : begin 

write (outf ile, t :7 :2, tab) ; 

for j :■ 1 to numTotEquations do 

if (j in printSet) then write (outf ile, u [j ]: 12 : sp [ j ], tab) ; 
writeln (outf ile) ; 
end; 
otherwise; 
end; (case) 

end; mm mm , 

{ Total function defining variables in terms of one another.™ ) 
function f(j: integer; t: extended; u: uArray) : extended; 
begin 


case j of 
xd 
x 

yd 

y 

wt 
th 


f 
f 
f 
f 
f 
f 

otherwise 


- thr*cos (u [th] /conv) *g (u [y] ) /u [wt ] ; 

- u[xd] ; 

- thr*sin (u [th] /conv) *g(u[y] ) /u[wt] - g(u[y]); 

- u [yd] ; 

* wtd; 

thd; 


end; 


end; 


begin 


Textbook (nil) ; 
tab ; " chr (09) ; 

for j :■ 1 to numTotEquations do 
sp [ j ] : - 2 ; 
sp [press] :■ 8; 
sptdens] :■ 8; 


open (outf ile, ' l/AE448prog.dat ); 
rewrite (outf ile) ; 
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m numlntEquations; t Number of integration equations } 

n trunc( (tend-tbegin) /dt) ; { Number of steps } 

{ Set names of all variables.- } 

printset [xd,x,yd,y, wt, vtot,th, gamma, thrust, press, dens] ; 

{ Set design variables— } 

wtO 2769714.29; { lbf } 

DryWtBooster : ~ 170000.00; { lbf } 

DryWtOrbiter 260000.00; { lbf } 

Time [Separation] :■ 148.00; { seconds } 

Time [Shutdown] :« 394.00; { seconds } 

Time [StartPitchl] 20.00; { seconds } 

Time [StopPitchl] 140.00; ( seconds } 

Time [ StartPitch2 ] 160.00; { seconds } 

Time [StopPitch2] :■ Time [StartPitch2] + 38.00; { seconds } 

[ Set initial conditions of all variables- } 
init Variables ; 

[ Echo initial conditions of all variables in the printing set... > 
print (3) ; 

{ Beginning of time-step loop— ) 
for i 1 to n do begin 

{ Computation of other, non time-dependent variables- } 

if (t m Time [Separation] ) then u[wt] u[wt] - DryWtBooster; 
thrMax thrustAvail (u [y] ) ; 

if (t < Time [Separation] ) then thr ;■ 10*thrMax else 
if (t < Time [Shutdown] ) then thr 4*thrMax else 

if (t < tend) then thr 0*thrMax; 

u[press] :« pressure (u[y] ) ; 
ujdens] density (u [y] ) ; 

u [thrust] ;■ thr; 

if (t < Time (Separation] ) then wtd :■ 10*wtd0 else 

if (t < Time [Shutdown] ) then wtd 4*wtd0 else 

if (t < tend) then wtd 0; 

thd pitchRate (t) ; 

if (u[th] > 360.00) then u[th] u[th] - 360.00; 
if (u[thj < -360.00) then u[th] :« u(th] + 360.00; 

( Runge-Kutta 4th-Order integration loop of all variables- ) 

for j :■ 1 to m do 

k (1, j] dt*f ( j,t,u) ; 
for j :■ 1 to n do 

uu( j] : “ u[ j] + k[l,j]/2; 
for j 1 to m do 

k[2,j] dt * f ( j , t+dt / 2 , uu ) ; 
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for j 1 to m do 

uu[ j] u[j] + k[2,j]/2; 

for j :» 1 to in do 

k[3, j] dt*f ( j,t+dt/2,uu) ; 

for j 1 to in do 

uu[ j] : m u(j] + k 1 3 , j ] ; 

for j 1 to m do 

k[4,j] dt*f ( j,t+dt,uu) ; 

for i 1 to m do 

u [ j] u[j] + {k[l, j]+2*k[2, j]+2*k[3,Dl+k[4 # Dl)/6; 

t :■ tbegin + i*dt; 

u[vtot] :■ sqrt (sqr (u [xd] ) +sqr (u [yd] ) ) ; 

if (u[xd] > 0) then u [gamma] :■ arctan(u[yd] /u[xd] ) *conv 
u[gamma] :■ 90.00; 

{ Printing parameters™ ] 

if (i mod cint - 0) and (printset <> []) then print(3); 


end; 

close (outfile) ; 


end. 


else 
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-162.16 
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-160.08 
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83.60 


83.01 
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80.58 


79.96 


79.34 
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8.556786 


8.329764 
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0.002063 


0.002037 
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0.001493 
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Figure 3.0.2: Plot of X-Position vs. Time 
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Figure 3.0.5: Plot of Weight vs. Time 
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Figure 3.0.6: Plot of Inclination vs. Time 
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Figure 3.0.7: Plot of Total Velocity vs. Time 
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Figure 3.0.9: Plot of Thrust vs. Time 
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Figure 3.0.10: Plot of Thrust/Weight Ratio vs. Time 
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Figure 3.0.1 1 : Plot of Y-Position vs. X-Position 
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4.0 Orbital Analysis 

The NASA ground rules for the Shuttle II specified two orbits into which the orbiter must be 
capable of launching: a 270 nautical mile orbit due east of Kennedy Space Center, and a 150 
nautical mile polar orbit To verify the feasibility of the proposed ASTS configuration, it 
was necessary to determine the types and magnitudes of the orbital maneuvers that would 
be necessary for orbit inserton. As a starting point, the output data for the final position and 
velocity from the launch simulation program was used. This data places the orbiter at an 
initial altitude of 78.1 nautical miles with a velocity of 25,962 feet per second. 


4.1 The Insertion Orbit 

The initial position and velocity of the orbiter place the orbiter originally at the perigee of 
an elliptical orbit. The semi major axis of this orbit, ai, was calculated using the Vis-Viva 

equation: . ^ ^ 

v‘" 

where ri is the sum of the altitude of the orbit and the radius of the earth (20925672.57 feet). 

The eccentricity of this first orbit , and the altitude and velocity at the apogee of this orbit 
were calculated using simple geometric relationships. The pertinent data on the first orbit 
is as follows: 


ei = .02473 
ai = 21,943,244 feet 
nap = 22,485,871 feet 
Vlap = 24,709 feet/second 


4.2 The 270 Nautical Mile Orbit 

In order to provide a minimum energy transfer between the insertion orbit and the 270 nmi 
orbit, an elliptical transfer orbit was used (see figure 4.2). 

The 270 nmi orbit has a radius of a2 = 22,567,272 ft. The semi major axis of the transfer 
ellipse, at, has a value of (ai +a 2)/2 = 22,255,258 ft. 


The velocities at periapsis and apoapsis of the transfer ellipse, as well as the velocity for 
the circular orbit were calculated using the Vis- Viva equation again. 


^ (t, ' a T ) 



7 


The following values were calculated for the velocities and the velocity changes needed for 
insertion into a 270 nmi orbit: 


Viap 24,709 ft/sec 

V t i 24,890 ft/sec 

dVi 181 ft/sec 


V t 2 24,799 ft/sec 

V 270 24,975 ft/sec 

dV2 176 ft/sec 

The total impulsice change in velocity for this orbital maneuver is 351 ft/sec. 


4.3 The 150 Nautical Mile Orbit 

The path of the elliptical orbit transfer orbit used for this maneuver is shown in figure 4.3. 

The velocity and radius of the 150 nmi orbit are found to be 

r 150 21837672 ft 

Viso ’ 25389 ft/sec 

Following the sane procedure as before in order to obtain the velocity changes necessary 
yields the following information: 

Vi a p 24709 ft/sec 

Vti 24678 ft/sec 

dVi -31 ft/sec 

Vt2 25419 ft/sec 

Vi 5 o 25389 ft/sec 

dV2 -30 ft/sec 


The total impulsive change in velocity for this orbital transfer is 
velocity changes would have to be made in a direction opposite to 
the shuttle. 


-61 ft/sec. These two 
the direction of travel 


of 
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Figure 4.3.1 : 150 Nautical Mile Polar Orbital Transfer 
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5.0 Stability and Control Analysis 

Stability and Control Analysis for the Advanced Space Transportation System was 
performed using geometric Shuttle data and empirical and graphical methods given by 
Etkin, Roskam and Smetana. The geometric parameters used to perform the analysis 
included physical data for the wing, the body, and the vertical tail. 

In order to determine stability and control derivatives, the geometric data for the Shuttle II 
must be used to determine aerodynamic parameters such as aspect ratio, wing taper ratio, 
and the mean aerodynamic chord. 

'A' is the aspect ratio and is defined to be equal to b^S; b is the wing span and S is the wing 
planform area. 'S' can be easily determined using the following formula: S=(c r +ct)(b/2), 
where c r is the wing root chord and ct is the wing tip chord. The ratio of the wing tip chord to 
the wing root chord is in itself and important physical parameter known simply as the 
taper ratio, X. The value of the mean aerodynamic chord (the wing chord which passes 
through the centroid of either half -wing) can be found by integrating the equation for the 
square of the wing chord: c 2 = {c r - [(c r -c t )/(b/2)]y) 2 over the span of half a wing and 
multiplying by 2/S. All of these aerodynamic parameters have been calculated and are 
found in the Executi ve Summary in Tables 3.2.1 and 3.3.1. 

5.1 Determination of Clcc 

CLa is an extremely important stability derivative describing the change in lift of an 
aircraft due to variation in angle of attack. 

Roskam gives the following formula for determining CLa : 

CLa = KwbCLaw 

where K w b =1 - .25(d/b) 2 + 0.25(d/b) 

and CLaw = 2jtA/(2+[(A 2 B 2 /K 2 )(l+tan 2 Ac/2/B 2 )+4] 1/2 } 

In these equations: d = body diameter at root chord = 30 ft. 

b = wing span = 127.7 ft. 

A = wing aspect ratio = 2.8 

P = compressibility parameter = (1 - M 2 )^ 2 (for subsonic flight) 

K = ratio of actual wing section left curve slope to 2n 
Ac /2 = sweep of the wing half-chord line = 34° 


Since these equations are useful only for Mach numbers less than the critical Mach 
number, K is assumed to be approximately 1. Hence, the following values for CLa are 
obtained: 


M range P? 

0.0 - 0.3 1.0 

0.3 - 0.6 0.775 

0.6-Mcr 0.415 


(/rad) 

2.97 

3.09 

3.33 


It is apparent that these values alone for the lift-curve slope (although correct) cannot be 
used to represent Cl<x for the ASTS. For this reason, methods of estimating CLa higher 
Mach numbers and varying Reynolds numbers (from Dynamics of Flight, by Etkin) were 
employed. 
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Etkin presents a graphical method by which the lift-curve slope of an aircraft may be 
obtained as a function of its geometrical parameters. Before incorporating Etkin s met , 
values for t/c and X (defined in Figure 5.1.1) were assumed to be 0.06 and 5 , respectively. 


Using Figure 5.1.1 the following section lift curve slopes were obtained: 


Bn 

lailflZlailoT 

iailQj/radl 

100 

0.87 

5.74 

10 7 

0.93 

6.14 

10 8 

0.95 

6.27 


It follows that these values for the actual lift curve slope are good approximations for 
transition both at the leading edge and at c/2. 

These values of (ai)o were then used in conjunction with Figure 5.1.2 to estimate the 
overall aircraft lift curve slope. Also, Etkin defines K as [(ai)om/2*]P and Ap as the 
compressibility sweep parameter equal to tan* 1 (tanAc/ 4 /p). 


Table 5.1.1. CLa of ASTS for M < Mcr 


M Range 

P 

Ab (°) 

Rn 

K 

pA/k 

pCLa/K </°) 

CLa(/rad) 

0.0 - 0.3 

1 

42 

io<J 

0.914 

3.06 

0.051 

2.67 

0.0 - 0.3 

1 

42 

i 0 r 

0.977 

2.81 

0.050 

2.80 

0.0-0.3 

1 

42 

lO 8 ^ 

0.998 

2.805 

0.049 

2.80 


0.88 

46 

10 s 

0.804 

3.06 

0.049 

2.57 


0.88 

46 

10 7 

0.860 

2.87 

0.048 

2.69 


0.88 

46 

10 8 

0.878 

2.81 

0.048 

2.74 

0.6 - Mcr 

0.644 

| 54 

10* 

0.589 

3.06 

0.0451 

2.36 

0.6 - Mcr 

0.644 

54 

10' 

0.629 

2.87 

0.043 

2.41 

0.6 - Mcr 

0.644 

54 

10 8 

0.643 

2.80 

0.043 

2.46 


Note that these values of Cl« are slightly higher than those estimated using the method 
outlined by Roskam. 

Etkin presents a system for predicting Cl« for Mach numbers greater than the critical 
Mach n um ber using the graphical method shown in Figure 5.1.3 which is a function o 
M, X, and Ad2- Results using this figure are listed in Table 5.1.2. 
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Table 5.1.2. CLa of ASTS for M > Mo- 


iMranae ! 

AtanAc /2 

A(M 2 -l) iy * 

CLa (/rad) 

ekbiH 

1.89 

1.64 

4.14 

EEBim 

1.89 

2.97 

3.42 

[m 

1.89 

r 4.04 

2.72 i 


1.89 

5.02 

1 2.21 


L89] 

f 5.97 

P““ 1.82 


1.89 

6.88 

1.62 



Prediction of CLa for Mach numbers pater than approximately 2.8 was not possible due to 
the lack of a proper estimation technique. 


5.2 Determination of the Neutral Point, hn 


The overall aircraft neutral point, hn, i. defined aa the point where the pitching moment is 
invariant with angle of attack. It ie further defined geometncally by Figure 5.2.1. 

Roskam gives the following relationship for determining the neutral point . 


hn = hnw + AhnB (for tail-less configurations) 


where hnw = ki[(XacVc r )-k2] 

and AhnB - - (dM/da)/(qScCLaw) 

where dM/da = (q/36.5) S Wf 2 (Xi)de/dctAXi 


Note that Roskam defines both the neutral point, hn, and the aerodynamic center, Xac/c, as 
the same quantity. 

r no effect on the 
point is reduced to 

the following form: 


Using the relationships for AhnB, it was found that the body had little c 
overall aircraft neutral point Therefore, determination of the neutral 


hn = ki[(Xac'/c r )-k 2 ] 

Xac'/cr ki and k 2 are found using Figures 5.2.2, 5.2.3, and 5.2.4, respectively Using 
ftfese figures together with known geometric data and the relationship for hn, the following 
results were obtained. 


Table 52.1. Neutral Point of the ASTS for M < Mcr 


M range 

P 


bi 

k 2 

hn 

0.0 - 0.3 

1.0 


1.43 

0.39 

0.27 

0.3 - 0.6 


0.59 

L43 

0.39 

0.29 

0.6 - Mcr 

0.64 

0.60l 

1.43 

0.39 

0.30 
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These values for hn show that as the Mach is increased, the neutral point begins to move 
gradually aft. This finding agrees with most aerodynamic data on similar aircratt. 

It is important to determine the characteristics of a flight vehicle in as many of its flight 
regimes as possible. Therefore, methods given by Etkin were used to find values of hn for 
Mach numbers greater than the critical Mach number, 

Etkin gives a graphical method shown in Figure 5.2.5 for determining hn as a function of 
A, MA and A d2- Using these variables and Figure 5.2.5, the following data was obtained. 


Table 5A2. Neutral Point of the ASTSfbrO< M<2£ 


M range 


■HBSEIS&Syi 

hn 

0.0 - 0.3 



0.28 


1.9 

2.5 

0.29 

lo.6 - Mcr 

1.9 

1.8 



1.9 

1.6 


1 1.3 - 1.6 

r i-9 

30 



1.9 

4.0 

o.4a 

HUM 

1-9 

5.0l 

0.48 

E OE23H 

1.9 

6.0 

r ~ 0.49 


1.9 

S 6.9 

0.49 


Again, these numbers show that (as expected) the neutral point starts out relatively close to 
the quarter-chord for low subsonic Mach numbers and moves aft approaching the halt- 
chord at higher supersonic Mach numbers. 


5.3 Determination of Cmo 

C mo is an extremely important stability derivative which dictates whether or not a 
configuration has positive or negative static stability- Cmo is determined through the 
following equation: 


Cmo = CLct (b - hn) 


where h represents the non-dimensional location of the aircraft center of gravity. For ® 
configuration to be statically stable, Cmo must be negative; this requires that the quantity 
(h - hn) must be negative. Xcg was determined by considering the location and weight of 
the components of the ASTS. The distance between the center of gravity and the leading 
edge of the mean aerodynamic chord (non-dimensionalized by dividing by the mean 
aerodynamic chord) is defined as h. Upon completion of these calculations, it is found that 

h = 0.136. 


Hence, Cmct = Cj/x (0.136 - hn) 


5.4 Clo of Vertical Tail 
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The lift curve slope of the vertical tail is determined by the geometric characteristics of the 
vertical tail and is influenced by the change in flow charactenstics due to the presence of 
the wing and body combination. Figure 5.4.1 shows the relationship between (ap)B/ai for a 
vertical tail on a body of circular cross section. Using a value of 1 for D/h, (aF)B/ai is 
approximately 2.9. 

Values ofcEafor the vertical tail are determined using the same procedure as before for the 
wing. When- this value is obtained it is multiplied by the factor (aF)B/ai to correct for its 
locationmthe flow field relative the the wing and body. It is assumed that the vertical tail 
has approximately the same values of t/c and t as the wing so that both the wing and 
vertical tail have the same values for (ai) 0 . Using the geometric parameters for the 
vertical tail and Figures 5.1.2 and 5.1.3 the following is found. 


Table 5.4.1. Corrected Clo values for the ASTS Vertical Tail 


M range 

Bn 

Cta 

Cl a corrected (rad* 1 ) 

0.0 - 0.3 

10 6 

2.93 

8.20 

0.0 - 0.3 

10 7 

3.08 

8.60 


10** 

3.14 

8.80 

lo.3-0.6 

”” —W 

2.62 

7.30 

EEBiXifli 

id 1 

2.74 

7.70 


10^ 

2.75 

| 7.70 

0.6 - Mcr 

10** 

2.36 

6.60 

0.6 - Mcr 

id 1 

P " 2.46 

6.90 

0.6 - Mcr 

W 

2.46 

— — 6 9Q 


— 

' 4.58 

” 12.80 

mmsm 

... 

4.35 

12.201 


Determination of Cut corrected for the vertical tail beyond a Mach number of 
approximately 1.6 was not possible. 


5.5 Determination of C| P 

Cio is the change in rolling moment due to a variation in rolling velocity. Ci p is often 
referred to as the rolling moment damping derivative. Etkin provides a graphical method 
by which Ci p can be calculated based on the following parameters: P, KA*,and Ae (which 
is defined as tan* 1 [(l/pltanAc/^). Etkin’s method is presented here in Figure 5.5.1. Table 
5.5.1 contains the values of Ci p found using Etkin’s method. 
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Table 5.5.1. Cip values ofASTS for M < Mcr 


M range 

“Pi 

Ae (°) 

Rn 

K 

Cm (/rad) 

0.0 - 0.3 


42 

1^ 

0.914 


0.0 - 0.3 

“ 1.0 

^ 42 

10^ 



0.0 - 0.3 

1.0 

42 

10 8 

0.998 

-0.230 

0.3 - 0.6 

0.88 

46 

10^ 

0.804 

-0.215 

0.3 - 0.6 

0.88 

46 

W 

0.860 

-0.225 

[•EffEfll 

0.88 

46 

lo 8 ! 

0.878 

-0.229 


0.644 

54 

10^ 

0.589 

-0.201 

0.6 - Mcr 

0.644 

54 

10 V 

0.629 

-0.210 

0.6 - Mcr 

0.644 

54 

10 8 

0.643 

-0.215 


C]p does not vary greatly with Mach number, and is always negative. 


5.6 Determination of Cip 

Cm represents the change in rolling moment due to a change in sideslip angle, p. Qp is 
primarily governed by the wing dihedral angle, the body and the vertical tail For this 
reason, Cip is often referred to as the effective dihedral derivative. Etkin present a 
graphical method of determining Cip based on the effects of the wing, the body and the 
vertical tail. This method utilizes Figure 5.6.1. Using this figure, it was found that 
Cipw = *.32 Cl. 

The effect of the body on C]p is found using the equation: 


(ACip)B = 1-2 (A) 1/2 Zw (h+wVb 2 

where Zw = vertical distance of c/4 point below fuselage centerline = 11.0 ft. 
h = average fuselage height at wing root = 28.0 ft. 
w = average fuselage width at wing root = 30.0 ft. 


Using this information: (ACip)B = 0.08. 

The effect of the vertical tail on Cip is found using the equation. 


(ACip)p = -ap (1 - do/dp)(Vp/V)2 Sp Zf/(S b) 

where Zp = vertical distance between vertical tail AC and aircraft CG = 41.0 ft. 
Sp = Vertical Tail planform area = 603 ft 2 . 

Using this equation under the assumption that da/dp is negligible and V F /V is 
approximately equal to 1: (ACip)F = -0.033 ap 


Combining components: Cip= -0.32 Cl + 0.08 - 0.033ap 

This equation holds for .til Mach numbers below the critical Mach number. 
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5.7 Determination of Cir 


C lr represents the change in rolling moment due to yawing velocity, or simply rolling 
moment due to yaw. Etkin divides Ci r into two components: Ci r due to the wing an ir 
to the vertical tail. Using Figure 5.7.1, (C] r ) B = 0-23 C L • 


due 


The contribution of the vertical tail on Ci r is shown through the following formula. 


(Cjr )f = aF Sf Zf [z +(lf/b)]/(S b) 

where If = distance between vertical tail AC and the aircraft CG = 48.0 ft. 


Hence, (Q r )F = 0-079 aF 


Combining components: C] r = 0.23 Cl + 0.079 ap 

This equation holds for all Mach numbers below the critical Mach number. 


5.8 Determination of C n p 


C n B is defined to represent the change in yawing moment due to sideslip angle, p. C n p is a 
very important stability derivative since it alone determines whether or not an aircraft has 
positive weathercock stability; for positive weathercock stability, C n p must be positive. 


Figure 5.8.1 shows the variation of C n (j/C L 2 as a function of wing aspect ratio, A, and wing 
quarter-chord sweep angle, A. 


Using Figure 5.8.1, C n p = 0.07 Cl 2 

5.9 Determination of C n p 

C np is a cross derivative defined as the change in yawing moment due to rolling velocity. 
Etkin presents the following equation for determining Cn p : 

C np = [(AC np )l/ClJ Cl + KACnpWDOla] (Cdo)<x 

where (AC np )l/CL and (AC np )2/(CD0)a are found from Figure 5.9.1 

From Figure 5.9.1, (AC np )l/CL = -0.17 and (AC^Cdo* = 19 > *> that C np may be written 

3 S 

Cnp = -0.17 Cl + 19 (Cdo)cl 

5.10 Determination of Cnr 

Cnr is defined to be the change in yawing moment due to yawing velocity and is often 
referred to as the damping in yaw derivative. Etkin presents the following equation for 
determining C n ri 

Cnr = KAC nr )i/CL 2 ]CL 2 + [(AC nr )2/CDo]CDO 
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where [(ACnrh/CL 2 ] and [(AC nr )2/CD0] are found from Figure 5.10.1 

Figure 5.10.1 reveals that [(ACnrh/CL 2 ] = -0-02 and [(ACnr^CDo] = -0.48, so that C nr may be 
written as 

Cm. = -0.02 Cl 2 - 0.48 Cdo 

5.1 1 Determination of CMq 

C mq is defined as the variation in pitching moment due to pitch velocity or simply as the 
damping in pitch. Smetana gives the following relationship for Cmq : 

C mq * - 2 x* lx' I Clo/c 2 

where x' = distance from CG to wing quarter-chord (positive for forward CG) = -6.97 ft. 

After substituting the proper geometric and stability parameters, C m q can be written as 
Cmq = 0.0372 CLa 

5.12 Determination of Ciq 

CLq is defined as the change in lift coefficient due to pitch velocity. Smetana gives the 
following relationship for CLq • 

Cj^q — 2 x Cjja/c 

Upon substitution of known geometric parameters, CLq is written as 
CLq = -0.2727 Clo 

5.13 Determination of Cd<x 

CDa is defined as the change in drag coefficient due to a change in the angle of attack. 
Smetana gives the following relationship for CDa: 

CDa = (dCDO/doO + [2 Ci/frA)] Clo 

Since (dCoo/da) is usually small, CDa may be written simply as [2 Cl/(*A)] CLa- 

Therefore, CDa may be expressed as 
Cdo = 0.2274 Cl CLa- 

5.14 Determination of Cyp 

Cyft is defined as the change in side force, y, due to a change in sideslip angle, P- Smetana 
breaks Cyp up into 3 components: that of the wing, body, and vertical tail. 

(Cyp)w = Cl 2 (6 tan A sin A)/[MA + 4 cosA)] 

Upon substitution, (Cyp) w = 0.2232 Cl 2 


60 



(C y p)B = -ki (Clo)b (Sb/S) 

where Sb = body reference area = 3733 ft^ 

ki = interference factor determined from Figure 5.14.1 


Figure 5.14.1 gives ki as a function of Z*/(d/2) where Zw is the distance from 1 body 
centerline to quarter-chord point of exposed wing root chord (positive for quarter-chord 
point below body centerline) and d is the max body height at wing/body intersection. 


Zw=ll-0ft. 
d = 25.0 ft 


From Figure 5.14.1, ki = 1.44 

Smetana gives the following values of (Clo)B ; 

(Cl<i) circular body = 0.0525 rad' 1 
(Cl<x) rectangular body = 0.1253 rad' 1 

Since the ASTS configuration is a combination of these two shapes, let (Cl<x)b - 0.08 rad 

Now that all values needed are known, (C y p)fi = -0.074 

(Cyp)F = -k (Clo)f (1 + da/d[3) Sp/S 

k is found from Figure 5.14.2 and is a function of (W(2ri)]. 

by = vertical tail span = 29.0 ft. 

ri = fuselage diameter in tail area = 30 ft. 

From Figure 5.14.2, k = 0.75 

If it is assumed that do/dB is negligible then (Cyp)F may be expressed as 
(Cyp)F = -0.078 aF 

Combining all components; Cyp = 0.2232 Cl^ - 0.074 - 0.078 aF 


5.15 Determination of C yp 

Cyp represents the change in side force due to rolling velocity. Smetana expresses C yp as 
the sum of the contributions due to the wing and the vertical tail. 

(Cyp ) w = Cl ([(A + cosA) tanMA + 4 cosA)] + 1/A) 

Substitution into (Cyp )w yields (Cyp )w = 0.91 Cl 

(Cyp )f = -2 aF Sf Z F /(S b) 

All of these parameters have been defined previously so that (Cyp )F = -0.067 aF 
Combining these two; Cyp = 0.91 Cl - 0.067 aF 
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5.16 Determination of Cyr 

Cyr is defined as the change in side force due to yawing velocity. Smetana expresses Cyr 
asthe sum of the contributions due to the wing and the vertical tail. 


(Cyr)w = 0.143 Cl * 0.05 


(Cyr)p = -2 If (Cyr)F/b = -0.752 (Cyr)F 
So that, Cyr = 0.143 Cl - 0.05 - 0.752 (Cyr>F 


5.17 Dynamic Stability Matrices 

Once the proper stability derivatives have been determined, it is possible to use them to 
investigate the dynamic stability of the configuration. Presented here are the longitudinal 
and lateral stability matrices corresponding to the non-dimensional, controls fixed, 
linearized equations of motion for a rigid aircraft. 


Longitudin al Eouations: 

Dx = Ax where x = (u a q 0) T 

A(l,l) » (Cxu + 2 Clo taneoV2n 
A(l,2) = Cxo/2p. 

A(l,3) = 0.0 

A(l,4) » -CL0/2p 

A(2,l) * (Czu * 2 CloV(2p - Czo) 

A(2,2) = Czo/(2|X - Czo) 

A(2,3) = (2 p + CzqV(2p - Czo) 

A(2,4) = -Clo tan©o/(2p - Cz a) 

A(3,l) = (C M u + Cmo A(2,l))/iB 
A(3,2) = (CMa + C mo A(2,2)ViB 
A(3,3) = (Cftfq + Cmo A(2,3))/iB 
A(3,4) = C ma AC2,4ViB 
A(4,l) = 0.0 
A(4,2) = 0.0 
A(4^) = 1.0 
A(4,4) = 0.0 


where p = m/(p S c/2) 

iB = iyy/tp s (c/2) 3 ] 

Cxa = Clo * CDa 
Cza = *Clo * Cdo 
Czu = * (M 2 /l-M 2 ) C lo 
CZq = -CLq 
CDm = CDu/M 

Cxu = -2 (Cdo + Clo tanQo) - M CDm 
Czo = 0.0 
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T .lateral Equations; 

Dy = By where y = (B p r <J>)^ 

B(l,l) = Cyp/2p 
B(l,2) = Cyp/2n 
B(l,3) = -(2 \i - CyrV2|i 
B(l,4) = Clo^IA 
B(2,l) = Gc Cip + iE C n pVA 
B(2,2) = (ic Ci p + i£ C np )/A 
B(2,3) = Gc Clr + iE CnrVA 
B(2,4) = 0.0 

B(3,l) = Ga Cnp + iE QpVA 

B(3,2) = Ga C np + iE Cl p )/A 

B(3,3) = GACnr + iEClryA 

B(3,4) = 0.0 

B(4,l) = 0.0 

B(4,2) = 1.0 

B(4,3) = tan6 0 

B(4,4) = 0.0 

where p = m/(p S h/2) 

A = iAiC*iE 2 
iA = Ixx/[p S (b/2) 3 ] 
ic = Izz/Ip S GV2) 3 ] 
ig = Ixz/[p S (h/2) 3 ] 

Once the matrices A and B are found, their Eigen Values can be determined. These Eigen 
values reflect whether or not the system is stable. For a system to asymptotically stable, all 
of its Eigen values must have negative real parts. 

The matrices A and B were calculated using a computer program written in FORTRAN. 
The program and results of the investigation of the dynamic stability for landing 
(assuming 0 O = 10° and U 0 = 300 ft/sec and L = W) are listed on pages 66 and 71. 


5.18 Determination of C|§ a 

Ci8a describes the variation in rolling moment with change in aileron deflection; often, 
this s tab ility derivative is referred to as the aileron power. Smetana gives the following 
equation for Cis®: 

Ab 

Ci5a = [2CLaX/(Sb)]\ cydy 

-'A 

where c is the equation for the wing chord as a function of wing span, y. 


c = c r - [(c r -c t )/(h/2)] y = 73.125 - 0.863 y 

and a and b represent the spanwise distance for the body centerline to the inner and outer 
aileron points, respectively. 



a = 27.0 ft b = 58.5 ft 

x is a correction factor which is a function of the aileron chord to wing chord ratio, t is 
determined from Figure 5.18.1. Since the ratio of aileron chord to wing chord is not 
constant over the span of the aileron, the average value was used to determine x. 

(c a /c w )avg = 0.33 t=0.57 

Upon integration and substitution of the proper values, 

C] 5ct = 0.071 CLa. 

5.19 Determination of Cisc 

Cj5r is the variation in rolling moment coefficient with change in rudder deflection. 
Smetana gives the following formula for estimating C]gr • 

Qfir = ap t Sp Zp/(S b) 

Here, t is a function of the rudder area to vertical tail area ratio. 

Sr = 157.5/10) Sf = 603 ft2 

Hence, S^F = 0.261 so that from Figure 5.19.1, x = 0.48. 

Substitution into the equation for C]gr results in 
CiSr = 0.0158 ap. 

5.20 Determination of C n sr 

Cn5r is the variation in yawing moment coefficient with a change in rudder deflection and 
is often referred to as the rudder power. Smetana gives the following relationship for C n 5r 
C n Sr = -*F * Sf lf/(S b) 

As with CiSr, the value of x is determined from Figure 5.19.1 and has a value of 0.48 so that 
Cn5r may be expressed as 

Cnfir = -0.0187 ap 
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AE448 CG Locations/1 



A 

B 

C 

D 

1 

Component 

X-Locat ion 

Waiaht 

Y-Moment 

2 

Winer 

116.60 

28.300 

3,299,780 

3 

Vertical Tail 

140.00 

3,200 

448,000 

4 

Forward Body 

28.90 

27,000 

780,300 

5 

Mid Body 

80.00 

23,600 

1,888,000 

6 

Aft Body 

128.00 

13,200 

1, 689, 600 

7 

Main Enaines 

137.00 


3, 123, 600 

8 

OMS 

138.00 


r 510,600 

9 

RCS 

138.00 

1,600 

220, 800 

m 

Control Systems 

48.00 

4,000 

192,000 

Vfl 

Avionics 

22.00 

3,100 

68,200 

m 

Personnel 

35.00 

1,600 

56, 000 

m 

Stores 

1 42.00 

600 

25,200 

m 

Encrine Accessories 

125.00 

6,200 

775,000 

m 

Hvdroaen Tanks 

64.00 

24,000 

1,536,000 

m 

Oxvaen Tank 

121.00 

11,000 

1,331,000 

eh 

TPS 

96.30 

6,100 

587,430 


Totals 


180,000 

16,531,510 

m 

Centroids 

91.84 
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f' He '!< ;'< CONTROI S FIXED LONGITUDINAL.. AMD I.AIRRhI- 
C *** DYNAMIC STABILITY AND CONTROL PROGRAM 
C *#* FOR THE SECOND GENERATION SPACE 
C #** TRANSPORTATION SYSTEM. 

C ### WRITTEN BY i PHILIP BENEFIELD & DEL JOHNSON 


C 


C 


C 


DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 


PRECISI ON G . S » W » US » CBAR . UO . RHO . THET AO 
PRECISION B . A » C . E . C L 0 . CDO . CMU . CM A 
PREC I S ION CMAD » CMQ . CL A » CLQ . CDA . CDU 
PRECISION CL.B . CLP » CLR . CNB . CNP . CNR 
PRECISION CYB . CYP . CYR . LLONG.MULGNG. IBLONG 
PREC I SION CXA . CZA . SOUND » M . CZl.1 . CZQ . CDM 
PRECISION CXU » CZAD . ALONG (4.4) »LLAT 
PRECISION IALAT . ICLAT » IELAT » MULAT 
PRECISION DEL . ALAT (4*4) 


OPEN(UNIT=9»FILE= '5TAB1IN.DAT' » ST ATUS= ' OLD ' ) 

• OPEN < UN IT* 10. FILE* ' STAB 1 OUT .DAT" .STATUS*' UNKNOWN ' ) 


READ ( 9 . & ) G . S . W . US . CBAR . UO » RHO » THET AO 

READ(9»#> B.A.C.E. CLO . CDO 

READ ( 9 » # ) CMU. CMA. CMAD »CMQ 

READ ( 9 . *) CLA . CLQ . CDA » CDU 

"READ ( 9 » * ) CLB . CLP . CLR 

READ ( 9 » * ) CNB » CNP . CNR . CYB . CYP . CYR 


00 

0 

* 

c 


iw 
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FORMAT ( IX » A40 . ' = / »F16.6) 

FORMAT (////////. 1 OX . ' GEOMETRICAL * INERT I AL 
' CHARACTERISTICS: './/> 


X EQUILIBRIUM 


y 


WRITE < 10 . 20) 
WRITE ( 10 » 300 ) 
WRITE ( 10 . 300 ) 
WRITE ( 10 . 300 ) 
WRITE< 10.300 ) 
WRITE ( 1 0 . 300 ) 
WRITE < 10 . 300 ) 
WRITE ( 10.300) 


'GRAVITATIONAL ACCELERATION ( FT/SEC~2 ) ' t G 
'WING AREA (FT~2)'»S 
'WEIGHT (LB)'.U 
'WING SPAN ( FT ) ' » WS 

'MEAN AERODYNAMIC CHORD (FT) '.CBAR 

' AIRSPEED ( FT/SEC ) ' » UO 

'AIR DENSITY ( SLUGS/FT”3 ) ' »RHO 


WRITE (10. 300 ) 

WR ITE < 10.300) 
WRITE ( 10.300) 
WRITE < 10.300) 
WRITE ( 10.300) 
WRITE ( 1 0 . 300 ) 
WRITE ( 10 . 300 ) 


WRITE ( 10.30) 

WRITE (10. 300 ) 

WRITE ( 10.300) 

WRITE (10.300) 

WRITE (10.300) 

WRITE ( 10.300) 

WRITE ( 10.300) 

WRITE (10.300) 

I.IPTTP ( 1 n . 300 ) 'CDU 


'INITIAL THETA (RAD) ' .THET AO 


' IVY 

(SLUG-FT-2) ' .B 

' I XX 

( SLUG-FT~2 ) ' .A 

' IZZ 

(SLUG-FT"2) ' »C 

' IXZ 

(SLUG-FT"2) ' .E 

'CLO 

( ) ' .CLO 

'CDO 

( ) ' . CDO 

x , •- non-dimensional 

' CMU 

( ) ' . CMU 

' CMA 

( > ' » CMA 

'CMAD 

i ( ) ' . CMAD 

'CMC 

( ) ' »CMG 

' CLA 

( ) ' .CLA 

'CLQ 

( ) ' .CLQ 

'CDA 

( ) ' .CDA 


. rni i 


DERIVATIVES t ' . 
qc’KV.WM- 

of POOR 


page IS 

QUALITY 




a nn 
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FOR M A T < / / / / p 10 X ? ' L A T E R A L 


D I RECTI ON At. 


WRITE (10? 40) 
WRITE ( 10 ? 300 ) 
WRITE ( 1 0 ? 300 ) 
WRITEC 10*300) 
WRITE ( 10 ? 300 ) 
WRITE ( 10 ? 300 ) 
WRITE ( 10 ? 300 ) 
WRITE < 10 ? 300 ) 
WRITE ( 10 ? 300 ) 
WRITE ( 10 v 300 ) 


' CLB (/RAD) ' »CLB 
'CLP < /RAD)'? CLP 
' CLR </RAD)'?CLR 
' CNB < /RAD)'? CNB 
'CNF ( /RAD ) ' f CNP 
'CNR < /RAD) 'r CNR 
'CYB </RAD)'?CYB 
'CYP (/RAD)'?CYP 
'CYR ( /RAD ) ' ?CYR 


derivatives: 


/) 


5(c^C5{t CALCULATIONS FOR LONGITUDINAL MOTION 


C 


LL0NG=CBAR/2 ♦ 0 

MULONG=< W/G) /(RHOTSTLLONG) 

IBL0NG=B/(RH0#S*(LL0NG#*3) ) 

CXA=CLO-CDA 
CZA=-CLA-CDO 
S0UND=1 116*4 
M=UO/SOUND 

CZU=-( (M*M)/<1.0-M*M) >*CLO 
CZQ=-CLQ 

CXU=“ < 2^0T < CDO+CLOTT AN ( THETAO ) ) )-<M*CDM) 
CZAD=0 . 0 


:*LONG ( 1 ? 1 > = 
ALONG ( 1 » 2) = 
ALONG < 1 ? 3 ) 
ALONG ( 1 ? 4 ) 
ALONG ( 2 ? 1 ) 
ALONG < 2 p 2 ) = 
ALONG < 2 ? 3 ) 
ALONG ( 2 p 4) 
ALONG < 3 p 1 ) 
ALONG ( 3 p 2 ) 
ALONG ( 3 p 3) 
ALONG ( 3 ? 4 ) 
ALONG <4?1) 
ALONG ( 4? 2) 
ALONG ( 4v3) 
ALiJN(j v 4 ') 


( CXU + 2 . OTCLOTT AN < THETAO ) ) / ( 2 . 0*MUL0NG ) 
CXA/(2.0*MUL0NG) 

0.0 

-CL0/(2.0*MUL0NG) 

( CZU-2 . OTCLO ) / ( 2 .OTMULONG-CZAD ) 
CZA/(2.0*MUL0NG-CZAD) 

( 2 ♦ OTMULONG+CZG ) / < 2 ♦ OTMULONG-CZAD ) 

■ ( -CLO*TAN ( THETAO ) ) / < 2 . 0*MUL0NG-CZAD ) 
~(CMUTCMAD*AL0NG(2p 1) ) /IBLONG 
=( CMAtCMADTALGNG<2p 2) ) /IBLONG 
ss < CMQFCMADT ALONG ( 2 p 3 ) ) /IBLONG 
=CM ADTALONG ( 2 ? 4 ) /IBLONG 
= 0 . 0 
= 0.0 
= 1.0 
'= 0.0 


c *:!<* CALCULATIONS FOR LATERAL MOTION 

C 

LLAT =WS/2 . 0 

MULAT=< W/G) / ( RKOTSTLLAT ) 

IALAT -A/ ( RH0T.FaLLATr.K3 ) 

ICLAT =C/ ( RHCTSTLLATTT3 ) 

I.ELAT =E/ ( RH0TSTLLATTT3 ) 

DEL=< IALAT# ICLAT) -IELAT#*2 


C 


c 


10 

c 


15 

C 

100 


101 


200 

201 

c 


ALAT < 1 y 1 ) ~CYB/ ( 2 . OY-MULAT ) 

AL AT < 1 9 2 ) ~ C Y P / ( 2 * OfcMULAT ) 

ALAT ( 1 3 ) «- ( 2 * 0&MULAT-CYR >/(?-♦ 0#MULAT ) 
ALAT ( 1 » 4 ) -CL.O/ < 2 * 0*MULAT ) 

ALAT(29 1>=( ICLAT*CLB+IELAT*CNB>/DEL 
ALAT ( 2 9 2 ) = ( ICLAT*CLP+IELAT*CNP ) /DEL 
ALAT ( 2 9 3) - ( X C L A T * C L R + 1 E I.. A T * C N R ) /DEL 
ALAT (294) =0 ♦ 0 

ALAT (3 9 1 ) = ( IALAT#CNB+IELAT*CLB)/DEL 

ALAT ( 3 9 2 ) = ( I ALAT*CNP+IELAT*CLP ) /DEL 

ALAT ( 3 9 3 ) = ( I ALAT#CNR+IELAT*CLR ) /DEL 

ALAT (39 4) =0*0 

ALAT ( 4 y 1 ) =0 ♦ 0 

ALAT ( 4 9 2 ) =1 ♦ 0 

ALAT ( 4 9 3 ) =TAN ( THET AO ) 

- ALAT ( 474F-C/ 7tT~ 


WRITE ( 10 9 100 ) 

DO 10 J 1 9 4 

WRITE( 109200) ( ALONG ( J 9 K ) •> K=1.94) 

IF ( J ♦ NE ♦ 4 ) WRITE (109201) 

CONTINUE 


WRITE ( 10 9 101 ) 

DO 15 J=l94 

WRITE( 109200) ( ALAT ( J 9 K ) 9 K = l94) 

IF ( J ♦ NE . 4) WRITE (109201) 
CONTINUE 


FORMAT (////// 9 25X 9 ' THE STABILITY MATRIX CORRESPONDING TO ' 
THE NONDIMENSIONAL 9 CONTROLS FIXED 9 LINEARIZED t ' * / * 
LONGITUDINAL EQUATIONS OF MOTION FOR' 9/9 

'A RIGID AIRPLANE ' 9 //) , ... , 

■<////»25X» 'THE STABILITY MATRIX CORRESPONDING TO 9 / 


9/9 


20X9 
25X9 
35X 9 ' 
FORMAT 
20X9 
23X9 
35X 9 
FORMAT 
FORMAT 


THE NONDIMENSIONAL 9 CONTROLS FIXED 9 LINEARIZED 


LATERAL 
A RIGID 
(5X9'I 

( 3 X 9 


EQUATIONS 
AIRPLANE' 9 //) 

' 9 4 ( F12 ♦ o 9 6X ) 9 ' 
974X9 ' i ' ) 


OF MOTION FOR' 9/9 


CLOSE ( 9 ) 
CLOSE ( 10) 
STOP 
END 


\ 


68 



GEOMETRICAL * I N E R T I A L 


& EQUILIBRIUM 


CHARACTERISTICS 


a ♦ 


GRAVITATIONAL ACCELERATION (FT/SEC 2) 

WING AREA (FT~2) - 
WEIGHT (LB) = 
WING SPAN (FT) = 
MEAN AERODYNAMIC CHORD (FT) = 
AIRSPEED (FT/SEC) = 
AIR DENSITY (5LUGS/FT~3) - 
INITIAL THETA (RAD) * 
IYY (SLUG-FT~2) = 
IXX (SLUG-FT"2) = 
IZZ ( SLUG-FT"2 ) = 
IXZ (SLUG-FT"2) = 

CLO ( ) = 

CDO ( ) = 


32 ♦ 174000 
5818* 000000 
180000,000000 
127.692000 
51 . 120000 
300.000000 
0.002377 
0.174500 

12490393.000000 
1686164,000000 

13399266.000000 
795369.000000 

0.418000 

0.050000 


non-dimensional stability derivatives: 

CMU ( ) = 

CMA ( — ) = 

CM AD ( ) = 

CMQ ( ) = 

CLA ( ) = 

CLQ < ) = 

CDA ( ) = 

CDU ( ) = 


0.000000 

-0.403200 

0.000000 

0.104200 

2.800000 

-0.763500 

0.266000 

0.000000 


lateral-directional derivatives: 


CLB 

(/RAD) 

- 

CLP 

(/RAD) 

* 

CLR 

(/RAD) 

- 

CNB 

(/RAD) 


CNP 

(/RAD) 

= 

CNR 

(/RAD) 


CYB 

( /RAD ) 


CYP 

( /RAD) 

- 

CYR 

(/RAD) 

- 


— 0 . 338000 
-0 . 225000 
0.775500 
0.012000 
-0.070000 
-0 . 030000 
— 0 ♦ 706000 
-0,195800 
0.514200 



- 0*00315899 
-0 * 02743687 
0*00000000 
0.00000000 


- 0.05570903 

- 0.74066202 

- 0*04074188 

0.00000000 


TH 


n c- 


THE STABILITY MATRIX CORRESPONDING TO 
NONDIMENSIONAL s> CONTROLS FIXED * LINEARIZED s- 
LONGITUDINAL EQUATIONS OF MOTION FOR 
A RIGID AIRPLANE 


0.00480166 

- 0.09003117 

- 0.00745438 

0.00000000 


0,00000000 

1.02411888 

0.00192645 

1.00000000 


- 0 * 0 1320 4 5 / 
-0 . 00232787 
0 , 00000000 
0,00000000 


THE STABILITY MATRIX CORRESPONDING TO ^ 

THE NONDIMENSIONAL i CONTROLS FIXED » LINEAR! ZED » 
LATERAL EQUATIONS OF MOTION FOR 
A RIGID AIRPLANE 


- 0.01545018 

- 0.50320992 

- 0.04867205 

1.00000000 


- 0.95942552 

1,69903754 

0.09279544 

0.17629303 


0.03298353 

0.00000000 

0.00000000 

0.00000000 
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Figure 5.1.1. Lift curve slopes for two-dimensional incompressible flow. 



, per degree 



Figure 5.1.2.Lift-curve slop, for swept and tapered wings at speeds below critical Mach 




per radian 



Figure 5.1.3. Lift-curve slope for swept and tapered wings at supersonic speeds. 
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Figure 5.2.2. Aerodynamic Center Location of Lifting Surfaces. 



0 .2 A . A *8 l-O 


Figure 5.2.3. Aerodynamic Center Transformation Constant Ki 






Figure 5.2.5. Chordwise position of the mean aerodynamic center of swept and tapered 
wings at high speeds expressed as a fraction of the mean aerodynamic 


chord. 



<«F>» <*y>»r 
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A t , degrees 

Figure 5.5.1. C] p for straight-tapered wings. 
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*0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 


Figure 5.6.1. Qp for straight-tapered wings with no dihedral. 


82 







Figure 5.7.1. Charts for estimating C,, for subsonic incompressible flow. 


A. degrees A. degrees 



t 






Figure 5.9.1. Variation of (AC n ph/CL 

incompressible flow. 


and (AC np )2/(CD0)a with aspect ratio for subsonic 
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Aspect ratio, A 


Figure 5.10.1. Charts for estimating (ACnrVCL 2 and (ACnr)2/C D 0 


A, degrees 


K| 



d/2 


Figure 5.14.1. Values for wing 


fuselage interference factor. 
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Figure 5.14.2. Values for k as a function of the ratio of vertical tail span to fuselage 
diamter in the tail region. 
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Figure 5.18.1.Values for x as a function of aileron chord to wing chord ratio. 



Figure 5.19.1. Values for x hs a 


function of rudder area to vertical tail area ratio. 
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